function rest = sgrdsin(t,y,T,alphaK,s)
%function rest = sgrdsin(t,y,T,alphaK)


msvalue=[s.m1(t,T),s.m2(t,T),s.m3(t,T),s.m4(t,T),s.m5(t,T),...
    s.m6(t,T),s.m7(t,T),s.m8(t,T),s.m9(t,T),s.m10(t,T),...
    s.m11(t,T),s.m12(t,T),s.m13(t,T),s.m14(t,T),s.m15(t,T),...
    s.m16(t,T),s.m17(t,T),s.m18(t,T),s.m19(t,T),s.m20(t,T),...
    s.m21(t,T),s.m22(t,T),s.m23(t,T),s.m24(t,T),s.m25(t,T)];
rest(alphaK+1)=0;
rest(1)=0;
rest(2:alphaK+1)=s.sigma(t)*msvalue(1:alphaK);
rest=rest';